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The study of post-Einsteinian metric extensions of general relativity (GR), which preserve the 
metric interpretation of gravity while considering metrics which may differ from that predicted by 
GR, is pushed one step further. We give a complete description of radar ranging and Doppler 
tracking in terms of the time delay affecting an electromagnetic signal travelling between the Earth 
and a remote probe. Results of previous publications concerning the Pioneer anomaly are corrected 
and an annually modulated anomaly is predicted besides the secular anomaly. Their correlation is 
shown to play an important role when extracting reliable information from Pioneer observations. 

\Q • The formalism developed here provides a basis for a quantitative analysis of the Pioneer data, in 

order to assess whether extended metric theories can be the appropriate description of gravity in 

f^*) , the solar system. 
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O . I. INTRODUCTION 

O; 

Experimental tests of gravity show a good agreement with General Relativity (GR) at all scales ranging from 
laboratory to the size of the solar system 0, IS 0- However there exist a few anomalies which may be seen as 
challenging GR. Anomalies in the rotation curves of galaxies or in the relation between redshifts and luminosities can 
be accounted for by considering dark matter and dark energy but they can as well be thought of as consequences of 
modifications of GR at galactic or cosmological scales 0, El . 

The anomalous acceleration recorded on Pioneer 10/11 probes might point at some anomalous behaviour of gravity 
at a scale of the order of the size of the solar system 0, 0]- The observation of such an effect has stimulated a 
significant effort to find explanations in terms of systematic effects on board the spacecraft or in its environment but 
this effort has not met success up to now The Pioneer anomaly remains the subject of intensive investigation 
because of its potential implications in deep space navigation as well as fundamental physics |lfl 1 1 1 L Il2l ll."'| . New 
missions have been proposed [l4| and efforts have been made for recovering data associated with the whole duration 
O of Pioneer 10/11 missions and submitting them to new analysis [Tolll^ . 

The present paper follows up publications which have investigated whether or not metric extensions of GR had 
, the capability to account for the Pioneer anomaly while remaining compatible with other gravity tests performed in 
the solar system. Such 'post-Einsteinian' extensions preserve the very core of GR with gravity identified with the 
metric tensor <7 M „ and motions described by geodesies. In particular, the weak equivalence principle, one of the most 
accurately verified properties in physics, is preserved. However the extended metric may differ from its standard 
(GR) form so that observations may show deviations from standard expectations. An important point is that these 
extensions explore a broader family of metrics than in the usually considered PPN family , including in particular 
deviations in the outer solar system. 

These extensions have been introduced in the context of a linearized treatment of gravitation fields [HOil and then 
discussed with non linearity taken into account [l9| . We will show that the previous studies were only preliminary 
and that the more precise and detailed investigations presented in this paper change some of their conclusions. But 
the main result, namely that the post-Einsteinian extensions of GR show the capability to account for the Pioneer 
anomaly, will not be affected. Objections to this statement, contained in recent publications [23, 0], will be shown 
to miss their target. 

The theoretical motivations for extensions of GR are rooted in its long confrontation with Quantum Field Theory. 
Their discussion, presented in [18l Il9|. is not repeated in the present paper. Here we will focus our attention on the 
phenomcnological implications of these extensions, by testing the metric in the solar system through its confrontation 
with observations, particularly those associated with Doppler tracking of Pioneer 10/11 probes. These Pioneer data 
show an anomalous acceleration ap directed towards the Sun with a roughly constant amplitude over a large range 
of heliocentric distances 

a P ~ 0.8 run s~ 2 , 20 AU < r P < 70 AU (1) 

Note that the positive sign for an acceleration directed towards the Sun has been chosen to fit the convention of 8]. 
The numbers are given as indications which will allow us to discuss orders of magnitude later on. The symbol AU 
stands for the astronomical unit. 
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Besides this secular term, the recorded anomalous acceleration also shows diurnal and annual modulations 
As the secular one, these modulated anomalies could be the consequence of some not yet understood artefact. But 
the search for an artefact accounting for the secular anomaly is usually focused on systematic effects on board the 
spacecraft or in its local environment and it is clear that modulated anomalies can certainly not be explained in this 
manner, since nothing in the vicinity of Pioneer probes is expected to have diurnal or annual variations. This entails 
that secular and modulated anomalies can hardly be due to the same artefact. 

The main result of the present paper will be that modulated anomalies are a natural prediction of post-Einsteinian 
extensions of GR. As a matter of fact, the Doppler observable not only depends on the motion of the Pioneer probe 
but also on the perturbation of electromagnetic propagation along the up- and down-links. As the paths followed 
by these links are themselves modulated by motions of the stations, the anomalous Doppler acceleration is expected 
to contain diurnal and annual modulations. The diurnal and annual anomalies have to be considered as further 
observables of great interest to be confronted to theoretical expectations. As these observables can be correlated with 
the more frequently discussed secular anomaly, this opens new perspectives for testing the metric in the solar system, 
even if the systematics associated with modulated and secular anomalies are likely not correlated to each other. 

In the following, we will give a common description of the secular and modulated parts of the anomalous acceleration 
by introducing a representation of the Doppler tracking observables in terms of propagation time delays. The advantage 
of this representation will be to treat in a natural and consistent manner the influence of metric perturbations on 
probe motion on one hand, link propagation on the other hand. The benefit will appear clearly in the discussion 
of Doppler observables, deduced by differentiation of the so called radar ranges, that is to say time delays between 
emission and reception, as well as in the interpretation of observations. We have to stress at this point that Pioneer 
10/11 missions were not equipped with range measurement capabilities, which is quite unfortunate. This indeed leads 
to ambiguities in the determination of ranges and will be shown to play an important role in the interpretation of the 
anomalies. 

Basic definitions and relations between the various quantities will be written down in the context of the 'post- 
Einsteinian' extensions of GR in the next section (sect. 2). The delicate problem of taking into account motions of 
Earth and probe will then be addressed (sect. 3). Exact relations will be presented as well as analytical approximations 
accurate enough for the purpose of the present paper. We will use the fact that the deviation of the extensions from 
GR certainly remains small since most gravity tests are compatible with GR (sect. 4). Using these theoretical tools, 
we will study the Pioneer-like anomalies possibly arising in Doppler tracking of probes in the solar system (sect. 5). 
In order to discuss the relevant orders of magnitude, we will then make simplifying assumptions, considering the case 
of probes moving in the ecliptic plane and having nearly radial motions in the outer solar system. We will present 
theoretical expectations for the secular anomaly as well as for the modulated anomaly due to the motion of the Earth 
(sect. 6), taking into account the correlation arising between these two anomalies. We will then draw conclusions 
(sect. 7) from the results of this new analysis. 



II. BASIC DEFINITIONS AND RELATIONS 



As already discussed, the high accuracy of tests of the weak principle of equivalence allows us to focus our attention 
on metric extensions of GR. This does not mean that there are no violations of this principle but only that such 
violations are too small to account for the large Pioneer anomaly (of the order of one thousandth of the Newton 
acceleration at the place explored by Pioneer probes). 

We also disregard the effects of rotation and non sphericity of the Sun which have an influence in the inner solar 
system but hardly in the outer one. Hence, we consider a static and isotropic metric representing space-time around a 
punctual and stationary source. This assumption notably simplifies the description with metric fields only depending 
on two functions <?oo and g rr of a single variable, the radius r, 

ds 2 = g Q0 {r)c 2 dt 2 + g rr (r) (dr 2 + r 2 d6 2 + r 2 sin 2 9dip 2 ) 

<?oSV = deg^ = d^g^u = (2) 

The metric has been written with Eddington isotropic coordinates recommended by the IAU convention [2^, [2j| to 
represent coordinates in the solar system |24|. Radii are defined from the gravity source (r = at Sun center), 
colatitude angles 6 are defined with respect to the ecliptic plane (9 = 7f/2), and azimuth angles if describe rotation 
within the ecliptic plane. For simplicity, we consider the Earth center of motion to have a uniform circular motion at 
frequency fli = 2w yr~ x . We also disregard the problems associated with diurnal rotation of Earth and atmospheric 
perturbations. These simplifications lead to the drawback that diurnal modulations will not be modeled. 

For the sake of precision, the following remarks have to be made with respect to the IAU convention. As already 
stated, we have first disregarded the effects of rotation and non sphericity of the Sun, since they have a small influence 
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on the observables studied thereafter. The idea is that differences between the values calculated with the standard 
metric of GR or with the modified metric @, are only slightly affected by this simplification. The preliminary 
evaluation of anomalies performed in this manner will have to be confirmed by more complete calculations taking into 
account perturbations due to the structure and rotation of the Sun as well as to the presence of planets Q . 

Then, the IAU convention J^] explicitly refers to GR whereas we are here considering extensions of GR. We then 
have to face the implications for the definition of fundamental constants, primarily the velocity of light c. To this 
aim, we write the metric components as sums of standard GR expressions and small deviations 

= [SV-L + S 9i*» , \ s 9^\ < 1 (3) 

and we convene that the deviations vanish at the radius of Earth orbit on which or in the vicinity of which the most 
accurate experiments are performed 

%„(ri)=0 , ri = fAU~f50Gm (4) 

It thus remains to study the effects of variations with the radius r of the anomalous metric components 5 goo and Sg rr . 
The standard metric is the GR solution written with Eddington isotropic coordinates 

W st = (^ T f) 2 . M* = -(i + £) 4 ( 5 ) 

The constant k is related to the Schwartzschild radius (with Gn the Newton constant) 

K =G^„ L5km „ 10 -8 AU (6) 

& 

The dimensionless potential n/r is small in the solar system, with a value ~ on Earth and even smaller values 
at the large radii explored by Pioneer probes. 

The extensions of GR are often discussed within the PPN framework [l( where the metric |J5J is expanded in terms 
of the Newton potential n/r and Eddington parameters f3 and 7 inserted in front of the terms of the expansion (with 
P = 7 = 1 in GR) 

2 

.900 = 1 - 2^ + 2/3^ + ... , g rr = -l-2-i--... , PPN (7) 

The PPN metric can be considered as a particular 'post-Einstcinian' extension of GR with anomalies showing specific 
dependences on the radius 

2 

^oo-203-l)^ , 5{goo9rr)^ -2( 7 -l)~ , PPN (8) 

In this paper, we consider more general extensions of GR, which show significant deviations at long ranges (outer 
solar system) and not only short ones (inner solar system). With respect to the PPN metric JHJ, the more general 
extensions can be thought of as allowing for anomalies in the two sectors which may depend on scale. 

Einstein curvatures corresponding to the extended metric have been studied in a detailed manner in 19]. In 
contrast to the standard expressions which vanish everywhere except on the gravity source, the anomalous curvatures 
are generally non null in space outside the source. This is already true for PPN extensions and, again, the general case 
corresponds to a more general r— dependence. This dependence can also be described in terms of running coupling 
constants which replace the Newton constant while depending on scale. We do not repeat these calculations [13 
but recall that the natural metric extension of GR involves two running coupling constants which correspond to the 
sectors of traceless and traced tensors [2^ ■ 

From the point of view of phenomenology, the two sectors are as well represented by the two functions 8goo (?*) and 
5 (goo9rr) (?")■ The first sector represents an anomaly of the Newton potential which has to remain small to preserve 
the good agreement between GR and gravity tests performed on planetary orbits Meanwhile, the second 

sector represents an extension of PPN phenomenology with a scale dependent parameter 7. It opens an additional 
phcnomenological freedom with respect to the mere modification of the Newton potential and this freedom opens the 
possibility to accomodate a Pioneer- like anomaly besides other gravity tests 0, 0, . 

Recent publications force us to be more specific on the relation between the Pioneer anomaly and modifications 
of the Newton potential, i.e. anomalies in the first sector according to the terminology of the preceding paragraph. 
Interpreting the Pioneer anomaly in such a manner requires that Sgqo varies roughly as r at the large radii explored 
by Pioneer probes. If this dependence also holds at smaller radii |8|, or if the anomaly follows a simple Yukawa 
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law [4|, one deduces that it cannot have escaped detection in the more constraining tests performed with martian 
probes [27], H^, |2!| . Brownstein and Moffat have explored the possibility that the linear dependence holds at distances 
explored by Pioneer probes while being cut at the orbital radii of Mars j^lj ■ lorio and Giudice [2£j as well as Tangen 
\2l\ have in contrast argued that the cphcmeris of outer planets were accurate enough to discard the presence of 
the required linear dependence in the range of distances explored by the Pioneer probes. This argument has been 
contested by the authors of [2O] so that the conflict remains to be settled. 

Iorio and Giudice [2(| and Tangen [2l| have pushed their claim one step farther by restating their argument as an 
objection to the very possibility of accounting for the Pioneer anomaly in any viable metric theory of gravity. Even 
before entering the detailed developments to be presented in this paper, we can show that this claim is untenable just 
because it only considers metric anomalies in the first sector while disregarding those in the second sector. Later on 
in the paper, we will come back to the discussion of the compatibility of the metric modifications with observations 
performed in the solar system, often with an accuracy higher than that of Pioneer observations. This has to be done 
with care, accounting for the presence of the two sectors as well as for possible scale dependences. This question has 
already been discussed in 0, 0] for the cases of deflection experiments on electromagnetic sources passing behind 
the Sun [nisi mnim We will see below that it has a particularly critical character for the ranging experiments 
which involve directly the Shapiro time delay |36| . 



III. RADAR RANGING AND DOPPLER TRACKING OBSERVABLES 



In the present section, we introduce the time delay function, a two-point function the knowledge of which is 
equivalent to a characterization of the metric. We deduce from this function the radar ranging observables between 
the Earth and a probe in the solar system. We then analyze the case of Doppler tracking observables which are 
obtained from ranging ones through a time differentiation. 

Starting from the static and isotropic metric written in terms of Eddington isotropic coordinates, we define the 
time delay |3?1 | as the time taken in this coordinate system by a light-like signal to propagate from a spatial position 
xi to another one X2 

x a = r a (sm8 a cos tp a , smOaSmtfa, cos 8 a ) , a =1,2 (9) 

This defines a two-point function T which depends on the positions xi and X2 only through three real variables, which 
can be chosen as the two radii r-y and r 2 and the angle <f> between the two points as seen from the gravity source, 

cost/) = cos 6\ cos 82 + sin 0i sin 62 cos (y> 2 — <fix) (10) 

Stated differently, the time delay T is a function of the triangle built on the emitter, the receiver and the gravity 
source. 

The form of the time delay function for a static isotropic metric was obtained in |l8l Il9| by solving the Hamilton- 
Jacobi equation for a light ray 

c T(rx,r 2 ^) = / , = , ^ , Pdr/r , (11) 



goo v ' r-* y goo v ' r J 



These quantities are integrals over the light ray of integrands depending only on the conformally invariant ratio 
Soo/Srr- The parameter p, hereafter called the impact parameter, is an implicit function of the variables r , 1 ,7 , 2 ,</> 
determined by the second equation in (|1 lfl . This definition fixes the relative sign of p and 4>. 
One can also deduce from the time delay function 



cT = p4>+ T C(r)dr , C(r) = J -^(r) - £ (12) 
Jri V 500 r 

Calculating the partial derivatives of T, one indeed notices that cancellations appear in the angular derivative, due 
to the particular form (|ll|l of the time delay function, 

cd ri T=-C(n) , cd r2 T = C(r 2 ) , c9 T = p (13) 

Second order derivatives of the time delay function T may then be written 

[I "V,, j (»> 

cd <j ,d r2 T= =d r2 p , cd 2 T = d r C{r 2 ) + 
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Similar relations hold involving derivatives with respect to ri instead of r 2 . Note that non local expressions only come 
from angular derivatives, namely p and d^p. It also turns out that the angular second derivative of the time delay 
function is positive (for r 2 > ri). 

We now study the observables which are used in deep space navigation |3S| and make explicit their relation to 
the time delay function We begin with radar ranging observables obtained by timing radio signals exchanged 

between stations on Earth and the deep space probe. An up- link radio signal is emitted from Earth at (i^x^ - ) and 
received and sent back by the probe at {p2, X2) with the down- link radio signal received on Earth at (if , x|~ ) (positions 
identified with the center of motion of Earth, itself assumed to run uniformly on a circular orbit). Meanwhile the deep 
space probe is assumed to follow a geodesic trajectory in the outer solar system. As the gravity fields are very small 
there, this trajectory can be approximated with a good approximation as a Keplerian hyperbolic trajectory escaping 
the solar system ||. 

The positions of the emission and reception events are connected by light cones, that is also by the following 
relations between event times and the time delay function, 

i 2 -*r =T(xr,x 2 )=T_ , t+ -t 2 =T(x+ x 2 ) =T+ (15) 

For probes equipped with range measurement capabilities (which was not the case for Pioneer 10/11), the ranging 
observable may be defined as half the time elapsed on Earth from the emission time t± to the reception time 

M^ t l^ = ^± (16) 
2 2 V ' 

This quantity is not a proper gauge invariant quantity but it is directly related to such a quantity, the proper time 
sf — sj" elapsed on Earth between the same two events, through a mere multiplication factor determined by the 
potentials created by Sun on Earth and the velocity of Earth on its orbit. This multiplication factor is constant in 
the simple context studied in this paper, and it can therefore be omitted. The ranging time is available for observers 
on Earth even if they don't have access to the transponding time i 2 defined on board the deep space probe. However, 
the transponding time can be deduced from the solution of equations of motion, at soon as the metric is known with 
a sufficient accuracy. We will see below how to deal with this complication. 

For Pioneer 10/11 probes, the tracking technique was based on the measurement of the Doppler shift, a proper 
observable y defined from the ratio of cycle counting rates of reference clocks located at emission and reception stations 

. dst dtt 

V = ln-^ = ln-i; (17) 
as 1 at 1 

Note that y has its definition gauge invariant when written in terms of proper times, but can as well be written in 
terms of coordinate times on Earth for the same reason as in the preceding paragraph. The same information can be 
encoded in a Doppler velocity 

V dAt dtt -dtT , y , 1 + Vic 

— = = -4: = tanh - , y = In '— (18) 

c dt dtf + dt± 2 y l-V/c y ' 

The difference between the observables V and y appears only at third order in the velocities and can be considered 
as a small term. This is due to the choice of the median observer time t to parametrize the data 

t + +t - 

t = 1 - 1 , dt~l = dt + dAt , dt 1 =dt- dAt (19) 

Should another time be used in its place, we would obtain second order corrections in the relation between V and y. 

The following remarks are worth being kept in mind when using tracking data for obtaining knowledge on the 
motion of the deep space probe. Clearly the Doppler velocity is primarily correlated to the velocity of the deep space 
probe relatively to that of Earth, with relativistic as well as gravitational corrections fully accounted for in relation 
(|18fl . But the velocity of the probe at some specific time is not known with a sufficient accuracy, unless informations 
extracted from Doppler data are used. It follows that the Pioneer gravity test is more appropriately discussed in 
terms of the Doppler acceleration A, which is just the time derivative of the Doppler velocity [8[, 

This observable A gives a more direct access to the acceleration law of the probe to be compared with the theoretical 
expectation. Now the expected acceleration depends on the distance of the probe to the gravity source and the 
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latter also suffers an imprecise determination on probes for which ranging data are not available. These 'ambiguity' 
problems will be faced in the following, keeping in mind that relativistic and gravitational corrections can in principle 
be affected by the discussion. 

In order to compare tracking observations with theory, we need to make more explicit the relation between the 
ranging observable At and the time delay function T. To this aim, we note that, as the Earth has been assumed 
to have a circular motion, the time delay function T reduces to a function of two variables r 2 and (j). For emission 
and reception events connected by light cones l|15|) . these two variables may equivalently be replaced by the time 
parameters t\ and t 2 which parametrize the Earth and probe trajectories respectively, at least once the metric and 
the trajectories are known. The history of ranging observations, i.e. of series of time triplets (see eas ll5|) . may then be 
conveniently represented under the form of a relation At(t) with At and t defined by l|ltjfl and (|19fl . The transponding 
time t 2 thus appears as a function of the median time t. These implicit relations may for instance be solved by 
approaching in an iterative manner the solution of the following equations 

T_=T(x 1 (t-At) ) x 2 (t a )) , T + =T( Xl (t + Ai),x 2 (t 2 )) (21) 

h = t+ — t , At = 



2 ' 2 
Derivatives of the function T with respect to t x and t 2 are deduced from (|T3|) 



cd 1 T = c^ = pd 1 cf> , cd 2 T = c^^C{r 2 )r 2 + pd 2 <P (22) 
oti ot 2 



The Doppler velocity and the derivative of transponding time are given by similar relations 

c + V 1 + d 2 T + 1 + d 1 T_ 2 1 - 2 T_ 1 + d 2 T + 



(23) 



c-V 1 - diT + 1 - d 2 71 ' d t t 2 1 + <9iT_ 1 - d Y T ^ 

IV. RANGING AND DOPPLER ANOMALIES 

In order to bring the relations written in the preceding section to explicit formulas, we have to face rather complicated 
expressions, which can be dealt with in a numerical procedure but hardly in an analytical calculation. A simpler 
approach, extremely useful for a first discussion of the anomalies, is to use a first-order expansion of the observables 
in the metric perturbations |19j. The basic methods to be used in such an approach are presented in this section. 

Let us first discuss the case of the time delay (|11H with its standard form 

[cT] at {r u r 2 ^) = if(r)A (24) 

At first order in the metric perturbation, the post-Einsteinian time delay may be written as the sum of this standard 
form and of an anomaly (see and l|12|) ) 

c5T{r l ,r 2 , < f) = {cT-[dT] at }{r 1 ,r 2 ^) = f* 5w(r) ~ [gr . r]st (r)-. 



[ffoo] Bt V '[C(r)] s 



Sw ^ 5 { = -fc- - (25 ) 
I V foo J 2[gy r ] st 2 [ ff00 j st 

We again notice that this variation is determined by the perturbation of the conformally invariant ratio g rr /9oo- in 
order to compute the variation of the ranging time observable <|15[) induced by that of the time delay l|25|l , we now 
introduce notations for anomalies of the ranging time T and transponding time t 2 

SAt = At-[At] Bt , St 2 =t 2 -[t 2 ] st (26) 

and write a first order equation for these anomalies by linearizing (|21|) 

[d,r +]st + [ftr_ ]rt _^ + (l + [d2 r +]st - [d2 r^ ^ _ st + - st_ 



2 V 2 / 2 

[fltTJrt - [9xT_] st ^ sM [d 2 T + ] st + [d 2 T_] st 5t = 5T+ + ST- 



(27) 
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The symbols appearing in (|27|l are defined according to 

5T ± = <5T( Xl (*±[A*] st ),x 2 ([t 2 ] st )) (28) 
[d a T±] st = [d a T] st (xi(t± [At] 8t ),x 2 ([t 2 ] 8t )) , a = 1,2 

Solving these equations, one obtains the ranging and transponding time anomalies 

6 At = 



l-[V] st /c ST+ l + [V] st /c 5T„ 



2 l-[9iT + ] st 2 l + [axT_] st 

The implicit equations 1)21(1 and other expressions written up to now will be approximated in the next section so that 
they are more easily used for explicit discussions of Pioneer-like anomalies. The approximation will be based on the 
fact that the Earth and probe velocities are much smaller than light velocity. 

We now come to the fact that the anomalous time delay between emission and reception is not only affected by 
the perturbation of light propagation, but also by the perturbation of the probe trajectory. Precisely, the value of r 2 
to be used in preceding calculations is not the same in extended theory as it would be in standard theory. This is 
characterized by differences (5r 2 , S9 2 , S(p 2 ) of the coordinates of the probe. In order to discuss this point, we introduce 
a new notation for variations of any quantity / as sums of terms, the first one being calculated with endpoints fixed 
and the other ones associated with position differences, 

5f ee /(r 2 , 0) - [/] st ([r 2 ] st , [<j>]J = Sf + 5r 2 3 r [f] st + tyfy [/] st (30) 

The angular anomalies have been collected in 5<fi and no anomalies have been accounted for the position of Earth 
which is assumed to be known. The condition 10} of null metric anomalies at Earth ensures the consistency of this 
description with the conventions of metrology. Note that, in a first order expansion in the metric perturbation, all 
contributions to l|3l)[l except the first one may be calculated using the standard expression of /. For the time delay 
function in particular, equation l|30[) is read as 

ST = ST([r 2 ] st , + [C] 8t {[r 2 ]JSr 2 + [p} st 8<j> (31) 

The perturbation ST defined by (|3*T)l and evaluated at first order will be used in the following as a good approximation 
to the ranging time anomaly (|29|l . Its evaluation still requires the solution of the equation of motion of the probe. 

In particular, the position differences Sr 2 and S4> have to be deduced from the geodesic equations written separately 
for the standard metric [g] st and the modified one [<?] st + 5g. These geodesic equations have their usual form in a 
metric theory [3!| 



ds 

with T^ p the Christoffcl symbols and the relativistic velocities 



I*yV = (32) 



g Xp dx p 
r ^ = —i.d^v P + d v g^ p -d p g^) , u p = — , g^u^u" = 1 (33) 



For the computations to be performed in the next section, it is worth noticing that the geodesic equations (|32|l may 
also be written as conservation laws for energy and angular momentum. The latter is a vector, and the conservation of 
its direction just means that the motion takes place in an orbital plane containing the Sun. This plane is characterized 
by two angles, the longitude of the ascending node Q and the inclination of the orbit L, which also give the spatial 
direction of the conserved angular momentum along the unit vector (sin 1 sin f2, — sin 1 cos f2, cos l). 

In order to describe motion, we then introduce angular coordinates adapted to the orbital plane with 9 = ^ on 
the orbit, and (p measured in the orbital plane. The reduced energy e and angular momentum j are defined as the 
following conserved quantities 

_ cdt ._ 2 d<f> 
e = goo— , j = g rr r — (34) 
as ds 

These relations just give the velocity components u° and , with the third component u 6 vanishing and the fourth 
one u r given by velocity normalization l|33|) 



dr _ r Trr /TTr , 2 _}-((* j 



:2 



z- = u r = ir , (try = [ 9oo U-^Sf ~ e ( 35 ) 

ds goog rr \ \ g rr r z 
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U r is a function of the variable r giving the radial velocity u r and depending on the form of the metric and on the 
conserved quantities labelling the trajectory. Equations (|34l) can also be written in terms of functions U° and of 
the variable r 

cdt n e dip - j , 

— = U° = — , -f = ^ = ^— 36 
as goo as g rr r 2 

Standard trajectories are obtained by integrating these equations, with the metric components having their standard 
expressions. At first order in the metric perturbation around GR, the probe trajectory is then obtained as the sum 
of standard and anomalous contributions. The latter are expressed in terms of variations of the functions U . In 
particular, the radial velocity and acceleration show the following anomalies 

^(S) =HU r d r Ur) = ld?{[U r } 2 st }6r+ 1 -d r 5{(U r ) 2 } (37) 



ds 2 

For the explicit calculations to be performed in the next sections, we will write the solution for the distance variation 
as 

Sr 2 = Sr* + [U r (r 2 )) st f ^dr (38) 



[U r ] 



5U r _ S(g og r 



[U r ] st 2 [.900.9, 



TJst 



r 2 \g rr J r 2 [g rr \ st 



goo9rr (U r f 

The constant Sr„ represents the initial radius variation at r* between geodesies calculated for the extended and 
standard metrics. The angular variation is then written as 



6U* Sj Sg r 



(39) 



The constant Sip* represents an initial angular difference between the geodesies. 

When taken with the results of the preceding section, these equations provide us with an exact description of radar 
ranging and Doppler tracking, in the simplified context considered in the present paper and in a first order expansion 
in the metric deviation from its standard form. 



V. PIONEER LIKE ANOMALIES 



Implicit equations written in the preceding section have to be solved in an iterative manner, which is well adapted 
to a numerical procedure but not easily performed in an analytical work. In order to be able to present qualitative 
but explicit discussions of Pioneer-like anomalies, we now introduce approximated forms of these equations. We will 
also consider Pioneer-like probes with high excentricity orbits, so that it will be possible to neglect angular terms. For 
simplicity, we also consider that the probe moves in the ecliptic plane, i.e. that the inclination of the orbital plane 
vanishes (i = 0). 

The main argument pleading for these approximations is the fact that the Earth velocity is much smaller than light 
velocity Qiri/c ~ 10 . This entails that the change of time delay function due to motion of Earth during the time of 
flight of the signal is small. Furthermore, the parity of equations (|21|l and (|29|) leads to corrections induced by Earth 
motion appearing only at second order in Earth velocity ((fliri/c) 2 ~ 10 -8 ). Hence, the modifications of the ranging 
time anomalies (|29|) due to Earth motion may be ignored in a first discussion. Note that this statement applies to 
the anomalous part l|29l) of the ranging time and to the effect of motion during the ranging time only. It does neither 
hold for the standard ranging time which suffers an appreciable effect due to the Earth motion nor for the effect 
of Earth motion on a longer term. The latter effect is evaluated below and found to play a significant role in the 
interpretation of the anomaly. 
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As already discussed, the anomalous time delay is also affected by the perturbation of the probe trajectory. This 
is taken into account by using equation (|31[l which now gives the true anomaly of the ranging observable 

cSAt(t) ~ cST(t) = c5T(n, [r 2 ] st , [flj + [C] Bt ([r 2 ] st )5r 2 + [p] at 5<f> (40) 

The anomaly of Doppler velocity observable is then obtained from (|18|l and 125|) 

!v - [ ^{ra^ +18 ^ ,1 « fr2 } +lcwl -'* 2 

+ K.tf+w.pgj; (41> 

where we have introduced shorthand notations f 2 = t?-, <j> = ^ and Sr 2 = tj 2 , S<j> = The time derivative of 
the impact parameter appearing in l|41|) is given by 

(42) 

All contributions to the Doppler shift anomaly l|41|l are local except some contributions to the impact parameter 
anomaly 5p 



Sp . . I f r2 — [g rr ] st Sw dr Sr 2 



Note that explicit dependences of C on r 2 are omitted from now on. 
The anomaly of Doppler acceleration observable is obtained similarly 



(43) 



OA - [r 2 J st r [r 2 J st C r I -y— ^ r^p- 

+ J t (N st <5r 2 ) + S^r 2 + [C] st 5f 2 



Wst \P\ st I 7 " 2 ] st I [grr] st SlV \ 

r Wrt ^+JU_s P \ (44) 



with shorthand notations f 2 = and <5f 2 = ^r 1 - The first and second line in l|44|l contain all the terms which are 
not modulated by the annual motion of Earth. In particular, they contain the secular contribution to the anomalous 
Doppler acceleration which was calculated in [l8Lll9l |. Note the relative signs between terms in first and second lines of 
(|44|l . which correct an error made in [l^.[T^|. The third line contains modulated terms depending locally on anomalies 
of the probe trajectory and shown below to give negligeable contributions to Pioneer-like anomalies. The fourth line 
includes terms which depend on the non local anomaly dp. These important terms (see below) were ignored in |l8lll9j . 
They vary with the Earth motion around the Sun and determine annual modulations of the anomaly. 

The Pioneer anomaly has been recorded on deep space probes with high excentricity orbits and, therefore, nearly 
radial motions. For the sake of simplicity, we neglect from now on all the terms proportional to angular velocities 
or angular accelerations. As proper time relativistic corrections depend on the probe velocity squared r 2 , it is also 
possible to use the simplification ds ~ cdt. In this context, one deduces from (|35l) the radial acceleration read as the 
sum of standard and anomalous contributions 

c 2 <9 r ln [g 00 ] t [r 2 ] 2 
[r2]st ~ 2 [g rr ] " ~ ~~2^ [ffOOflrrU 

c 2 6(d r g 00 ) r „ , S(g 00 g rr ) 
or 2 ~ — -. : L r 2Jsti r~ 

* Woogrrlst Woogrrlst 

N s 2 t S{d r (g 00 g rr )) d r [gooffs] st 1 ( 

o f i f 1 L r 2j s t 0r 2 

* L.9003rrJ st Woogrrlst 
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We have disregarded the effect of planets as gravity sources, which is justified once again by the fact that we focus 
our attention on the anomalies (this effect has to be taken into account in the data analysis process ||). Equation 
(|35Jl is integrated to obtain the anomaly expressed on the position or velocity of the probe 



c Sgoo/2 - c 2 e st Se <%oo5rr) 

dr 2 - — p i — m tf r~ hit 



Sr 2 * [r 2 ] st Lr + f ^ = ^ ^ = W * S{9m9rr) dr) (46) 

{ Jri 2 [r] st [gooffrr] st J 

These relations are just particular cases of the more general results obtained in the end of the preceding section. They 
have been written in terms of a difference St of epoch, that is also the time of passage at the initial radius r*. For 
convenience, this initial radius has been pushed back to r\ . 

The Doppler acceleration anomaly (|44J) may then be expressed in terms of the radial velocity [ra]^ and acceleration 

c 2 [CL 6(d r g Q0 ) 



5A 



'2 



[>2 



[goog rr ] st 

[ c ] st %oo5rr) , S(g 00 g rr ) 



st [ [goog rr ] st 2 [ 300 ] s 2 t [C] st 
2 J [ C L | S(d r (g 00 grr)) d r [ 



f S(g ogrr)\ . 1 d [C] st (% 05rr) 



(9 r 



V 2 [5oo c W 2 N 8t * [goog 



1 d ^ st [C] ^ ^009rr] Bt \ C 2 6gc~ 




"•Jst 



N st * S [500ffrr] st J 2 [ 500 ] st 

K t [p] st N st r^ 2 og !Z )_ [ j 

Mst [ 5 oo] s 2 t [C] s 3 t I 2 

^ ( [r 2 ] st [0 r C] rt Sr 2 + [ P ] st ^ + j^hrSp \ (47) 



We now evaluate this expression by using orders of magnitude known for the different quantities. We first know that 
the metric components are close to 1 for 500, to -1 for g rr , the differences being of the order of the Newton potential 
n/r. The latter has a value ~ 10 -8 on Earth orbit and values 20 to 70 times smaller at the distances explored by 
Pioneer probes. The square of velocity divided by light velocity has the same value ~ 10 -8 for Earth, due to the 
virial theorem, and it is roughly 6 times smaller for Pioneer probes (velocity ~ 12 km/s ~ 0.4 times that of Earth). 
As we study deep space probes at large heliocentric distances r 2 r%, we also use the fact that the parameter [p] st 
is at least 20 times smaller than [r 2 ] st , so that terms scaling as [p]^ t / [r 2 ]l t are at least 400 times smaller than unity. 
For the same reason, terms proportional to angular velocity anomaly 5<j) are found to have a negligible effect. After 
these remarks, expression (|47|) is simplified to the following dominant contributions 

5 A ~ SA SCC + SA anu (48) 
SA sec ~ -— d r (6g 00 ) + [hU ( S ^oog rr ) _ gg \ _ c g2 ^ , ^ 



2 

d 



•^ann — ,, 

at 



st 



The term 8A SCC contains secular contributions proportional to metric anomalies in the first and second sectors as well 
as to the range ambiguity 

x r *i Sx P c2S 9oo ~ 2 ° 2 Nst 6e ~ lf] 2 s t S (9oo9rr) j \ ( . 

Sr 2 ^ [r 2 \ st <St- —3 S dr \ (49) 

{ Jrx 2[r] st J 

The last term SA ann is a modulated contribution which is proportional to the anomaly of the impact parameter 



h - - M st \pU + too ^ + (so) 
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At this point, it is worth emphasizing the differences between the expression (|48|l obtained for the anomalous accel- 
eration BA and the result previously obtained in 0, . 

The secular anomaly arising from the second sector replaces the result obtained in preliminary calculations 0, 0] 
which was spoiled by a calculation error. The previous result was linear in the gravity fields and proportional to the 
kinetic energy of the probe. The new expression appears at the second order in gravity fields, since it is proportional 
on one hand to the anomalous potential S(googrr) at the position of the probe and on another hand to the standard 
probe acceleration [^2] st , that is also the gradient of the standard Newton potential. It however remains of first order 
with respect to the metric anomaly, and this property has been used in the derivation. Note that the change of form 
of this term has no consequence on the comparison of the anomalies recorded on Pioneer 10 and 11 since the two 
probes had nearly equal velocities and anomalous accelerations . But it will affect several conclusions to be drawn 
in the next section. 

The previous results 0,^3 were also preliminary for the following reasons. First, the secular anomaly is corrected 
by a term proportional to the range ambiguity, because the position of the probe is not known directly (no range 
measurement capabilities on Pioneer probes), and this important fact was not discussed previously. The range 
ambiguity (|49J) contains contributions proportional to anomalies as well as trajectory mismodeling, i.e. modifications 
of the standard acceleration observable due to changes St and Se of the constants of motion. Then, the range 
ambiguity also affects the evaluation of the annually modulated anomaly, probably the most striking new feature of 
the expression (|48|l . Expectations for the annual and secular anomalies are thus correlated, a property which will 
turn out to be of uttermost importance in the next section. Finally, the comparison of these expectations may open 
a road to a genuine test of the 'post-Einsteinian' phenomenological framework, as discussed below. 



VI. DISCUSSION OF ORDERS OF MAGNITUDE 



Obviously, the interest of this new road depends in a critical manner on the orders of magnitude of the various 
terms, to be discussed in the present section. This discussion is heavily dependent on the stringent constraints put 
on possible metric anomalies by the gravity tests already performed in the solar svstem. In order to write down the 
relevant arguments, we introduce potentials <J>at and $p in the two sectors as in [l8L fl9| 

5oo ~ 1 + 2$jv , -9oo9rr ^ 1 + 2$ P (51) 

We also use the simplest form of the standard impact parameter [p] . which is modulated by Earth rotation 

[p] st ~ -r x sin - iconj)) (52) 

The time t CO nj corresponds to a conjonction (closest approach) of Earth and deep space probe. The simple expression 
is sufficient for the purpose of the present section. 
In this context, the secular and modulated anomalies in (|48|l are reduced to 



6A SCC ~ -c 2 d r 5<5> N {r 2 ) + + £$ P (r 2 ) + 26$ N (r 2 ) 

r 2 z I r 2 





Sr 2 


= ai — < 


r 2 


r 2 { 


r 2 
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G N M 


ri 2 


7*1 2 



dA 

8A ann ~ -_Acos(2Qx(t -t conj )) - ^-^sin(2f7 1 (t-i conj )) 



d\ ee — — = — —r — = flfri ~ 6 x 10 — ms~ (53) 
Ti z rx 

These expressions of the secular and modulated anomalies can be considered as the key predictions of the post- 
Einsteinian framework presented in this paper. 

The secular anomaly SA sec contains a first term describing the effect of a Newton law modification and a second one 
gathering the contributions of the range ambiguity and of the potentials. The annual anomaly 5A &nn is determined 
by an amplitude A which contains contributions of the range ambiguity and of the two anomalous potentials. This 
amplitude is not annually modulated but suffers a secular change during the probe's journey. Considering that A 
varies slowly over a year, the modulation in (|53|l is essentially at twice the orbital frequency. The observed behaviour 
has a richer structure but its discussion must take into account the following remarks. First the constants of 
motion, and therefore A, are modified at maneuvers which are occuring twice a year on the average. Then, the 
mismodeling contributions on angles ip 2 and 9 2 , not studied in detail here, have periods f2i and 2Qi , as soon as a non 
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null inclination angle t is accounted for. These remarks qualitatively explain why the observed annual modulation 
may be more complicated than the simple expression <|53[) . 

The maneuvers are not frequent for Pioneer 10/11 probes, which is one of the main causes of their excellent 
navigational accuracy They interrupt the free geodesic segments by changing the values of the constants of 
motion. This essentially amounts to a change of the local velocity with no change of the position if the maneuver 
can be modeled as nearly instantaneous. The detailed description of the maneuvers is one of the most delicate 
parts of the data analysis process and it certainly goes much farther than the purpose of the present paper. In the 
following paragraphs, we focus our attention on the main features which can be qualitatively expected rather than on 
quantitative results. To this aim, we can go along with the simplified expression (|53[1 . 

The hierarchy of magnitudes appearing in (|53l) . a 1 r 1 2 /r 2 2 in SA SCC and a 1 r 1 /r 2 in SA auni plays a central role in this 
discussion. This is the reason why we have used in l|53|) the acceleration cii of Earth on its orbit as the natural scale 
for measuring anomalous Dopplcr accelerations, except for the first term which is the direct effect of a Newton law 
modification. We briefly discuss this term now, before embarking on a more complete analysis of the terms induced 
by the range ambiguity and the second potential and putting into evidence the correlation between the secular and 
modulated anomalies. 

Should the Pioneer anomaly be explained by an anomaly in the first sector, a linear dependence of the potential 
8<S>n would be needed to reproduce the fact that the anomaly has a roughly constant value Q over a large range of 
heliocentric distances rp 

c 2 d r 6<5> N ~a P , 20 AU < r P < 70 AU (54) 

The simplest way to modelize the anomaly would thus correspond to a potential varying linearly with r and vanishing 
at Earth orbit to fit the convention (0J 

6$ N ~ (55) 
tp 

We have introduced a length £p characteristic of the Pioneer anomaly 

tp 1 = % ~ 0.8 x 10~ 26 m" 1 ~ 1.2 x 10~ 15 AU" 1 (56) 
cr 

Should this model effectively describe the metric in the vicinity of Earth and Mars, its effects could not have escaped 
detection in the very accurate tests performed with martian probes such as Viking |27J . Numbers sheding light on this 
point are given in ||| (see also where similar conclusions are obtained for a modification of the Newton potential 
having the form of a Yukawa potential). The effect of the perturbation Ij55(l on planets would produce a change of 
their orbital radius. The order of magnitude (|5fi|l would lead to range variations of ^50km and ^lOOkm respectively 
at smallest and largest distances. Meanwhile, the Viking data constrain these measurements to agree with standard 
expectations at a level of ^100m and ^150m respectively. These numbers are different enough to eliminate the simple 
model (15511 with the coefficient lp chosen to fit the Pioneer anomaly. Note that the effect of the Shapiro time delay in 
the range evaluation, which should in principle have been taken into account in the discussion, has here a negligible 
influence [IJ. 

As already discussed, these results do not prove that the Pioneer anomaly cannot be reproduced in a metric theory. 
First, there is the possibility that the linear dependence needed to reproduce (|54f) at distances explored by Pioneer 
probes is cut at the orbital radii of planets on which the strongest constraints are obtained 30] . As it was already 
discussed, it then remains to decide whether or not the ephemeris of the outer planets are accurate enough to forbid 
the presence of the linear dependence (|54(l in the range of distances explored by the Pioneer probes pfl l2l| . This 
point remains to be settled 30]. In any case, there is another possibility, namely that the Pioneer anomaly is induced 
by the second anomalous potential 8§p rather than the first one 5&n- We now consider these terms which are still 
here even if there is no anomaly at all in the first sector (8&n = 0). More thorough studies will have to be performed 
later on to study the correlated effects of anomalies in the two sectors. 

We now focus our attention on the secular anomaly 6A sec and the annual amplitude A which are determined in 
(|53|l by the range ambiguity and the second potential 



SA SCC ^ ai r 4\ 2 -^+ 6* P (r 2 )\ , A~ai— (— -r 2 f ^dr \ (57) 



ri 



With the same assumption S&n = 0, the range ambiguity is given by 



Y 2 ^7Y M ^r + / te dr _r>6*p d ) (58) 
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We have introduced the non relativistic reduced energy e which allows us to express the standard velocity [r] st as a 
function of distance, using Q35p. 



(e 2 - 1) 




(59) 



The secular anomaly and annual amplitude appear in (|5*7jl as different linear superpositions of the reduced range 
ambiguity 5r2/r2 and of the second potential 5&p. The term in front of the annual amplitude turns out to be larger 
than the term in front of the secular anomaly by a factor T^jTv. This could appear to be contradictory with the fact 
that the annual anomaly is only a fraction of the secular one [8| in the data, but this is not the case. As a matter of 
fact, the data analysis process described in 8] is based on the a priori assumption that there is no annually modulated 
anomaly in the physical signal of interest. As explained above, this assumption is valid for anomalies induced by 6$ m 
but not for anomalies induced by 5$>p. In the context of this assumption, the choice of the best trajectory fitting 
the data tends to produce a null or quasi null value for the annual anomaly. This corresponds to a choice of the 
trajectory leading to a perfect or nearly perfect compensation of the two contributions to A in l|57|) . in conformity 
with the assumed absence of annual anomalies. 

In the context of the present paper where annual anomaly of the Doppler acceleration is a natural expectation, the 
observations reported in Q have a different significance. They mean that the unknown epoch St characterizing the 
motion of the probe has been fixed so that the two contributions to A compensate each other at some well chosen 
radius T2 

\A\ «a^ , ^r 2 r^dr (60) 
T2 r 2 J ri r A 

As a consequence of already presented arguments, this compensation has to be effective within a fraction of the order 
of or even smaller than r\/r2- Now, this compensation cannot remain perfect over a long period of time. The first 
reason for that is the different dependence on T2 of the two terms to be compensated by each other. The second 
reason is due to the maneuvers which change the constants of motion and thus affect the compensation. A precise 
estimation of the annual anomaly thus requires a complete solution of the motion including a detailed description of 
the maneuvers. As this task is outside the scope of the present paper, we will not be able to conclude whether or 
not the annual modulations reported in Q are effectively accounted for in a quantitative manner by the effect of the 
second potential. 

Despite this deficiency, the description just given of the annual anomaly nevertheless leads to a quantitative esti- 
mation of the secular anomaly. Equation i|tj(Jfl indeed fixes the otherwise unknown range ambiguity Sr2/ V2, so that 
SA^ Pr can be rewritten 



2 



»'2 



f'l 



<L4 sec ~ am J { — I Xp(r)dr + X p(r2)} , 8<S> P {r) = X p{r)r z (61) 



A roughly constant anomaly is produced when xp is constant, i.e. when <5$p(r) is quadratic in r, in the range of 
Pioneer distances. Identifying the expression 5A SCC to the observed Pioneer anomaly a p fixes the value of the constant 

8A sec ~-a P -> xp^-T r V-- 4xl0_8AU_2 ( 62 ) 

6Klp 

This value is 3 times smaller from what would have been obtained without accounting for the range ambiguity. Note 
that xp can take different values outside the range 20-70 AU of Pioneer distances and that it is not even forced to be 
exactly constant in this range. In fact, we know that xp has to vanish at Earth radius in order to obey the convention 
0. We also show in the next paragraph that it may have to be smaller than 1)62(1 between Earth and Mars in order 
to be compatible with planetary observations. 

To this aim, we consider a simple model with the potential obtained as the sum of linear and quadratic terms 
vanishing at Earth orbit to fit (J3J 

M> P ~ Jr - n) 2 + Mr - n) 

Snip 

The quadratic coefficient has been fixed according to (|62|) . The further characteristic length fj,p has been introduced 
to represent the radial derivative of the metric anomaly at Earth orbit. This linear term has to be small enough to be 
dominated by the quadratic one at distances explored by Pioneer probes (fj,p -C rp). Now the metric perturbation 
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(|63|l has also an influence on the already discussed range measurements on martian probes. As a consequence of the 
Shapiro effect, a range variation is found with a value of the order of ^700m, in conflict with Viking data. This 
conflict may be cured by cutting off the simple dependence 1)63(1 at the orbital radii of Earth and Mars. It is easily 
checked out that this does not affect significantly the predictions made for the Pioneer probes which are at much 
larger distances. 



VII. CONCLUDING REMARKS 

As stated in the introduction, the present paper follows up publications [TtI ITsl IT^| devoted to the study of post- 
Einsteinian metric extensions of GR. A more complete theory of radar ranging and Doppler tracking has been given 
in terms of the time delay function, allowing us to discuss the annual anomaly besides the secular one. A mistake 
made in previous publications in the evaluation of the secular anomaly has been corrected and important new results 
have been obtained. In particular, the annual anomaly has been found to be correlated with the secular anomaly 
through terms arising from ambiguities in the position of the probe. This correlation, in principle available in the 
data, has to be scrutinized in order to extract reliable information from Pioneer observations. 

The following qualitative statements summarize the results of the present paper. As the secular anomaly, the 
annual anomaly is a natural consequence of the presence of a second potential. This has to be contrasted with the 
first potential which does not produce a significant effect along the links |4(j . Then, the annual anomaly produced 
by propagation along the up- and down-links can be compensated near an arbitrary point by an appropriate choice 
of the trajectory of the probe. In fact this compensation is an output of any best fit procedure based on the a 
priori assumption that there is no annual anomaly. However, the compensation cannot remain exact when the probe 
moves as the two compensating terms have different dependences on the heliocentric distance. It follows that the 
annual anomaly reappears, either after some free evolution or after the next maneuver. This qualitative behaviour 
is reminiscent of the observations of annual anomalies which were reported in [jj. This situation certainly pleads 
for pushing this study and comparing the theoretical expectations with Pioneer data. It is only after a quantitative 
comparison, taking into account all the details known to be important for data analysis [8j , that it will be possible to 
decide whether or not the post-Einsteinian phenomenological framework does fit the observations. 

It is clear after these remarks that some of the conclusions of our previous papers have to be amended: the 
secular anomaly turns out to be proportional to the standard acceleration and to the second potential. The corrected 
expression is quadratic, and no longer linear, in the gravity fields, with one contribution standard and the second one 
anomalous. Identifying the expectation with the observed Pioneer anomaly now points to a second potential with a 
quadratic dependence on the radius. This corresponds to a constant curvature with an unexpectedly large value in 
the outer solar system (eq. ljrJ2)l of the present paper) . This quadratic dependence may have to be cut off at distances 
exceeding the size of the solar system as well as in the inner solar system (in order to pass Shapiro tests on martian 
probes). Note that the discussion of the preceeding section was mainly focused on the change of the Shapiro delay, 
due to the anomaly on $p. The modification of the orbital radii, which could in principle play a role, has been ignored 
because it was expected to have a negligible influence. The evaluation of correlated effects of anomalies in the two 
sectors will be necessary in order to be able to reach definitive conclusions. 

These conclusions constitute motivations for new experiments in the solar system. Clearly, experiments with ranging 
capabilities will offer qualitatively better perspectives than Pioneer observations which were performed without such 
capabilities. Missions going to the borders of the solar system will either prove or disprove the existence of 
the anomaly at such long distances. Comparison with the theoretical expectation of the present paper will give an 
answer to the question whether such an anomaly may have a metric origin, with the metric possibly departing from 
the GR prediction. This idea could also be tested on a shorter time scale by adding specially designed instruments 
on planetary probes going to Mars, Jupiter, or Saturn, the reduction of the explored heliocentric distance being 
compensated by a potentially large improvement of the measurement accuracy. 
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